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Abstract 



O ■ Using Liouville space and superoperator formalism we consider pure stationary 

states of open and dissipative quantum systems. We discuss stationary states of 
open quantum systems, which coincide with stationary states of closed quantum 
systems. Open quantum systems with pure stationary states of linear oscillator are 
suggested. We consider stationary states for the Lindblad equation. We discuss 



CN \ bifurcations of pure stationary states for open quantum systems which are quantum 

^ ■ analogs of classical dynamical bifurcations. 

o . 

^ I Introduction 

O The open quantum systems are of strong theoretical interest. As a rule, any microscopic 
system is always embedded in some (macroscopic) environment and therefore it is never 
really closed. Frequently, the relevant environment is in principle unobservable or it is 



^ . unknown [1, 2]. This would render the theory of open quantum systems a fundamental 
^ ', generalization of quantum mechanics [3, 4]. 

qh' Classical open and dissipative systems can have regular or strange attractors [5, 6]. 

Regular attractors can be considered as a set of (stationary) states for closed classical 
systems correspondent to open systems. Quantization of evolution equations in phase 
space for dissipative and open classical systems was suggested in [7, 8] . This quantization 
procedure allows one to derive quantum analogs of open classical systems with regular 
attractors such as nonlinear oscillator [7, 9]. In the papers [7, 8, 9] were derived quantum 
analogs of dissipative systems with strange attractors such as Lorenz-like system, Rossler 
and Newton-Leipnik systems. It is interesting to consider quantum analogs for regular 
and strange attractors. The regular " quantum" attractors can be considered as stationary 
states of open quantum systems. The existence of stationary states for open quantum 
systems is an interesting fact [10]. 

In this paper we consider stationary pure states of some open quantum systems. These 
open systems look like closed quantum systems in the pure stationary states. We consider 
the quantum analog of dynamical bifurcations considered by J.M.T. Thompson and T.S. 
Lunn [12] for classical dynamical systems. In order to describe these systems, we consider 

^ E-mail: t arasov@t heory. sinp . msu . ru 



1 



Liouville-von Neumann equation for density matrix evolution such that this Liouville 
generator of the equation is a function of some Hamiltonian operator. Open quantum 
systems with pure stationary states of hnear harmonic oscillator are suggested. We derive 
stationary states for quantum Markovian master equation usually called the Lindblad 
equation. The suggested approach allows one to use theory of bifurcations for a wide 
class of quantum open systems. We consider the example of bifurcation of pure stationary 
states for open quantum systems. 



II Pure stationary state 

In the general case, the time evolution of the quantum state \pt) can be described by the 
Liouville-von Neumann equation 

||a) = a|a), (1) 

where A is a Liouville super operator on Liouville space, \p) is a density matrix operator 
as an element of Liouville space. For the concept of Liouville space and superoperators 
see the Appendix and [18]-[21]. For closed systems, Liouville superoperator has the form 

k^~{LH-RH) or A = L-, (2) 

where H — H{q,p) is a Hamilton operator. If the Liouville superoperator A cannot be 
represented in the form (2), then quantum system is called open, non-Hamiltonian or 
dissipative quantum system [20, 21, 13]. The stationary state is defined by the following 
condition 

A\pt) = 0. (3) 
For closed quantum systems (2), this condition has the simple form 

Lnlpt) = Rnlpt) or L]j\pt) = . (4) 

In the general case, we can consider the Liouville superoperator as a superoperator func- 
tion [7, 8, 20]: 

A = A(L3„L+) or A = A{Lx,Rx), 

where X is a set of linear operators. For example, X = {q,p, H} or X — {Hi, .., Hg}. In 
the paper we use the special form of the superoperator A such that 

s 

A = -^{Lh -Rh) + Y1 ^kNkiLn, Rh), 

k=l 

where Nk{LH, Rh) are some superoperator functions and is an arbitrary nonzero 
superoperator. 

It is known that a pure state is a stationary state of a closed quantum system (1) 
and (2), if the state is an eigenvector of the Liouville space for superoperators Lh 
and Rh- 

Lh\p^) = \p^)E, Rh\p^) = \p^)E. (5) 
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Equivalently, the state is an eigenvector of superoperators and such that 

Lh\p^) = \p^)E, Lfjlpii,) = \p^) -0 = 0. 
The energy variable E can be defined by 

E = {I\LhM = {I\RhM = {I\L+M. 
The superoperators Lh and Rh for hnear harmonic oscillator are 

It is known that pure stationary states = p%^ of linear harmonic oscillator (6) exists 
if the variable E is equal to 

K = ^^(2n + 1). (7) 

III Pure stationary states of open systems 

Let us consider the Liouville-von Neumann equation (1) for the open quantum system 
defined of the form 

^ \pt) = -^iLH- Rh) \pt) + Yl FkNkiLn, Rh) \pt). (8) 

fe=i 

Here F*^ is some superoperator and N^ILh, Rh), where k = l,...,s, are superoperator 
functions. 

Let is a pure stationary state of the closed quantum system defined by Hamilton 
operator H. If equations (5) are satisfied, then the state |p^,) is a stationary state of the 
closed system associated with the open system (8) and is defined by 

||Pt) = -^(^^-^^)lA)- (9) 

If the vector is an eigenvector of operators Lh and Rh, then the Liouville-von Neu- 
mann equation (8) for the pure state has the form 

J] Nk{E,E), 

k=l 

where the function Nk{E,E) are defined by 

Nk{E,E) = {I\Nk{LH,RH)\p^). 
If all functions N}.{E,E) are equal to zero 

Nk{E,E)=Q, (10) 
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then the stationary state of the closed quantum system (9) is the stationary state of 
the open quantum system (8). 

Note that functions Nf;{E, E) are eigenvalues and is the eigenvector of the super- 
operators Nk{LH, Rff), since 

Nk{LH,RH)M = MNk{E,E). 

Therefore stationary states of the open quantum system (8) are defined by zero eigenvalues 
of superoperators Nk{LH, Rr)- 



IV Open systems with oscillator stationary states 

In this section simple examples of open quantum systems (8) are considered. 

1) Let us consider the nonlinear oscillator with friction defined by the equation 

j^Pt = -l[H:Pt] - ^f3[q',p'pt + Ptp\ (11) 
where the operator H is the Hamilton operator of the nonlinear oscillator: 

Equation (11) can be rewritten in the form 

= L-M + 2mffL-,{±-(L;f + ^(Ltf - Ai,) |,.,, (12) 

where A — Q,^ — and the superoperator Lj^ is defined for the Hamilton operator H of 
the linear harmonic oscillator by (2) and (6). Equation (12) has the form (8), with 

N(Lh,Rh)^^{Lh + Rh)--^Li, F^2m(3L-,. 

In this case the function N{E, E) has the form 

N{E,E) = E-^. 

Let 7 = The open quantum system (11) has one stationary state of the linear 

harmonic oscillator with energy En — {hu;/2){2n + 1), if A = 2pfiu;{2n + 1), where n is 
an integer non-negative number. This stationary state is one of the stationary states of 
the linear harmonic oscillator with the mass m and frequency cu. In this case we can have 
the quantum analog of dynamical Hopf bifurcation [12, 14]. 

2) Let us consider the open quantum system described by the time evolution equation 

^\pt) = L-\pt) + L-cos(^^L+^ \pt), (13) 



where the superoperator Ljj is defined by formulas (2) and (6). Equation (13) has the 
form (8) if the superoperators F and N{Lh, Rh) are defined by 
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m=0 

(14) 

The function N{E,E) has the form 

iv(E,^)^c„.( ( ) . 

m=0 ^ ' 

The stationary state condition (10) has the solution 

£; = |(2n + l), 

where n is an integer number. If parameter Eq is equal to hw, then quantum system 
(13) and (14) has pure stationary states of the linear harmonic oscillator with the energy 
(7). As the result, stationary states of the open quantum system (13) coincide with 
pure stationary states of the linear harmonic oscillator. If the parameter Eq is equal to 
hjjj{2m + 1), then quantum system (13) and (14) have stationary states of the linear 
harmonic oscillator with n{k, m) — 2km -\- k + m and 

£^n(fe,m) = ^(2/c + l)(2m+l). 

3) Let us consider the superoperator function Nk{LH, Rh) in the form 

N.iURH) = ^Y.''knvl^{'2LlR^-Ll^"^-Rl+n. 

n,m 

and all superoperators are equal to L/. In this case, the Liouville-von Neumann 
equation (8) can be represented by the Lindblad equation [22, 23, 20]: 

7 • -1 

J^\Pt) = ~{Lh - RH)\Pt) + ^ - LyLy, - RylK) \Pt)- (15) 

3 

with linear operators defined by 

n m 

If |p^) is a pure stationary state (5), then all functions Nk{E,E) arc equal to zero and 
this state |p^,) is a stationary state of the open quantum system (15). 
If the Hamilton operator H is defined by 

In "^^^ 2 ^/ \ 

2m^ +^-? +^y^P + P(l): 

then we have some generalization of the quantum model for the Brownian motion of a 
harmonic oscillator considered in [11]. Note that in the model [11] operators Vk are linear 
Vfc = akP + bkq, but in our generalization (15) and (16) these operators are nonlinear. For 
example, we can use = + bkH"^. 
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V Dynamical bifurcations and catastrophes 

Let us consider a special case of open quantum systems (8) such that the vector function 

Nk{E,E) = {I\Nk{LH,RH)\p), 
be a potential function and the Hamilton operator H can be represented in the form 

s 

H = Y,Hk. 

k=l 

In this case we have a function V{E) called potential, such that the following conditions 
are satisfied: 

where Ek = {I\Lhi^\p) — {I\RHk\p)- If potential V{E) exists, then the stationary state 
condition (10) for the open quantum system (8) is defined by critical points of the potential 
V{E). If the system has one variable E, then the function N{E, E) is always a potential 
function. In general, the vector function Nk{E,E) is potential, if 

dNk{E,E) _ dNi{E,E) 
dEi dEk ■ 

Stationary states of the open quantum system (8) with the potential vector function 
Nk{E,E) is depend by critical points of the potential V{E). It allows one to use the 
theory of bifurcations and catastrophes for the parametric set of functions V{E). Note 
that a bifurcation in a vector space of variables E — {Ek\k — 1, s} is a bifurcation in 
the vector space of eigenvalues of the Hamilton operator Hk- 

For the polynomial superoperator function Nk{LH, Rh) we have 

N n 
n=0 m=0 

In general, m and n are multi- indices. The function Ni^{E, E) is a polynomial 

N 

Nk{E,E) = Y,a'^J:^E\ 



n=0 



where the coefficients are defined by 



m=0 

We can define the variables xi — Ei — ai {I — 1, s), such that functions Nk{E, E) 
Nk{x + a,x + a) have no the terms x'^~^. 

N N n ^ 

n=0 n=0 m=0 ' ^ ^' 
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If the coefficient of the term x?' ^ is equal to zero 



(fe)^^ (fc) , ik) ^ a^k) (k) (k) ^ Q 



then we have the following coefficients: 

(fc) 

- (ifc) ■ 

If we change parameters an \ then an open quantum system can have pure stationary 
states of the system. For example, the bifurcation with the birth of linear oscillator pure 
stationary state is a quantum analog of dynamical Hopf bifurcation [12, 14] for classical 
dynamical system. 

Let a vector space of energy variables be a one-dimensional space. If the function 
N{E, E) is equal to 

n-l 

N{E, E) = ±anE'' + J2 "i^^ ^ - 2> 
then the potential V{x) is defined by the following equation 

V{x) = ±x"+^ + J2 ^ ^ 2, 

and we have catastrophe of type A±n. 

If wc have s variables Ei, where / = 1,2, ...,s, then quantum analogs of elementary 
catastrophes A±n, D±n, E±q, E^ and Eg can be realized for open quantum systems. Let us 
write the full list of potentials V{x), which leads to elementary catastrophes (zero modal) 
defined by V{x) — Vo{x) + Q{x), where 

n— 1 



A±n- Vq{x) ^ ±x'l+^ + ^ajx{ n>2, 

71—3 n— 1 

D±n: Vo{x) = xlx2±xr' + Y,o^^A+ Yl 4"^""'^ 



j=l j=n-2 
2 5 

E±e ■ Vo{x) = {xl ± + ^ ajxi + ^ ajXix^'^, 

j=l j=3 

4 6 

E7 : Vo{x) = xl + xixl + ttjX-i^ + ajXix^^^, 

j=l j=5 

3 7 

Es : Vo{x) — xl + xl + ajx{ + ajXiX^"^. 

j=l j=4 

Here Q{x) is the nondegenerate quadratic form with variables X2, x^, Xg for A±n and 
parameters 2:3, Xg for other cases. 
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VI Fold catastrophe 



Let us consider the Liouville-von Neumann equation (8) for a nonlinear quantum oscillator 
with friction, where multiplication superoperators Lh and Rh are defined by (6) and 
superoperators F and N{Lh, Rh) are given by the following equations: 

F = -2L-L+ , N(Lh, Rh) = aoLj + aiL% + a2{L+)\ (17) 

In this case, the function N{E, E) is equal to N{E, E) = aQ + aiE + a2E'^. 

A pure stationary state of the linear harmonic oscillator is a stationary state of 
the open quantum system (17), if N{E, E) — 0. Let us define the new real variable x and 
parameter A by the following equation: 

4aott2 - Oil 

X — E + - — , A 



2^2' ■" 4q;2 

Then we have the stationary condition N{E, ii^) = in the form — \ = 0. If A < 0, 
then the open quantum system has no stationary states. If A > 0, then we have pure 
stationary states for a discrete set of parameter values A. If the parameters ai, ai and A 
satisfy the following conditions 

^hu;{n-{ h— ), X^hfu^ — , 



2a2 '2 2" 4 

where n and m are non-negative integer numbers, then the open quantum system (17) has 
two pure stationary state of the hnear harmonic oscillator. The energies of these states 
are equal to 

1 1 

En^hv{n+-), En+m^ fuv{n + m + -). 



VII Conclusion 

Open quantum systems can have pure stationary states. Stationary states of open quan- 
tum systems can coincide with pure stationary states of closed (Hamiltonian) systems. 
As an example, we suggest open quantum systems with pure stationary states of lin- 
ear oscillator. Note that using (8), it is easy to get open (dissipative) quantum systems 
with stationary states of hydrogen atom. For a special case of open systems, we can use 
usual bifurcation and catastrophe theory. It is easy to derive quantum analogs of classical 
dynamical bifurcations. 

Open quantum systems with two stationary states can be considered as qubit. It allows 
one to consider open n-qubit quantum system described by (8) as quantum computer with 
pure states. In general, we can consider open quantum systems as a quantum computer 
with mixed states [15, 16, 17]. A mixed state (operator of density matrix) of n two- 
level quantum systems (open or closed n-qubit system) is an element of 4"- dimensional 
operator Hilbert space (Liouville space). It allows one to use quantum computer model 
with four- valued logic [15, 16, 17]. The quantum gates of this model are real, completely 
positive, trace-preserving superoperators that act on mixed state [15, 17].. Bifurcations 
of pure quantum states can used for quantum gate control. 

This work was partially supported by the RFBR grant No. 02-02-16444. 
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Appendix 

Liouville space 

The space of linear operators acting on a Hilbert space ?i is a complex linear space Ti. 
We denote an element A oi T-ihy a, ket- vector \A). The inner product of two elements 
\A) and \B) of H is defined as {A\B) = Tr{A^B). The norm \\A\\ = ^/{A\Aj is the 
Hilbert-Schmidt norm of operator A. A new Hilbert space Ti with the inner product is 
called Liouville space attached to H or the associated Hilbert space, or Hilbert-Schmidt 
space. 

Let (I a; >} be an orthonormal basis of 7i: 

<.W>=S(.-.'), /rf*><x| = /. 
Then \x,x') — \ \x >< x'\) is an orthonormal basis of the Liouville space H: 

My,y') = Si.-my-y'), l^.l^.'MM = i. (is) 

For an arbitrary element \A) of 7i we have 
where 

{x,x'\A) = Tr{{\x >< x'l^A) = Tr{\x' >< x\A) =< x\A\x' >= A{x,x'), 

is a kernel of the operator A. An operator p of density matrix {Trp = 1, p^ = p, p > 0) 
can be considered as an element \p) of the Liouville space H. Using (19), we get 

\p) = J dx J dx'\x,x'){x,x'\p) , (20) 

where the trace is represented by 

Superoperators 

Operators that act on H, are called superoperators and we denote them, in general, by 
the hat. 

For an arbitrary superoperator A on H, which is defined by A\A) — \A{A)), we have 
{x,x'\A\A) = JdyJ dy'{x,x'\A\y,y'){y,y'\A) = J dy J dy'A{x,x',y,y')A{y,y'), 

where A{x, x', y, y') is a kernel of the superoperator A. 
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Let A be a linear operator in the Hilbert space Ti. We can define the multiphcation 
superoperators and Ra by the following equations 

La\B) = \AB) , Ra\B) = \BA). 

The superoperator kernels can be easy derived. For example, in the basis \x,x') we 
have 

{x,x'\La\B) = j dy J dy'{x,x'\LA\y,y'){y,y'\B) ^ j dy J dy'LA{x,x' ,y,y')B{y,y'). 
Using 

{x,x'\AB) =< x\AB\x' >^ J dy J dy' < x\A\y >< y\B\y' >< y'\x' >, 

we get kernel of the left multiplication superoperator 

{LA){x,x',y,y') =< x\A\y >< x'\y' >= A{x,y)6{x' - y'). 
Left superoperators are defined as Lie and Jordan multiplication by 

L^B = ^{AB - BA), L\B = -{AB + BA). 

The left superoperators and right superoperators are connected by 

= ~^A^ = ^A- 

An algebra of the superoperators is defined [8] by the following relations 

1) Lie relations 

^A B ^ ^A^B~ ^B^A- 

2) Jordan relations 

^(AoB)oC ^" ^B^C^A ~r ^A-^C^B — ^AoB^C + ^BoC^A + ^AoC^B-: 



^{AoB)oC + ^B^C^A + ^A^C^B ^C^AoB + ^B^AoC + ^A^BoCi 

f+t+ -Lf+f+ -L f+f+ f+ f+_Lf+ f+_Lf+ f + 

^C^AoB + -^S-^AoC + -^^-^SoC ^AoB^C + -^BoC"^^ -^^oC-^S- 

3) Mixed relations 

^A-B ^ ^A^B ~ ^B^A^ ^AoB ^ ^A^B + ^B^A^ 

where 

A- B = ^{AB - BA), AoB^]-{AB + BA). 
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